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Detection of SiO emission from a massive dense cold core 
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ABSTRACT 

We report the detection of the SiO (J = 2 — > 1) transition from the massive cold dense 
core G333. 125—0.562. The core remains undetected at wavelengths shorter than 70- 
/jm and has compact 1.2-mm dust continuum. The SiO emission is localised to the 
core. The observations are part of a continuing multi-molecular line survey of the giant 
molecular cloud G333. Other detected molecules in the core include 13 CO, C 18 0, CS, 
HCO+, HCN, HNC, CH 3 OH, N 2 H+, SO, HC 3 N, NH 3 , and some of their isotopes. In 
addition, from NH 3 (1,1) and (2,2) inversion lines, we obtain a temperature of 13 K. 
From fitting to the spectral energy distribution we obtain a colour temperature of 18 
K and a gas mass of 2 x 10 3 M Q . We have also detected a 22-GHz water maser in the 
core, together with methanol maser emission, suggesting the core will host massive 
star formation. We hypothesise that the SiO emission arises from shocks associated 
with an outflow in the cold core. 

Key words: stars: formation - ISM: molecules - ISM: jets and outflows. 
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1 INTRODUCTION 

1.1 Dense cold core G333. 125-0. 562 

The massive, dense, cold core G333. 125— 0.562 (a./2ooo = 
16' 1 21 m 35 s ,(5j2ooo = -50 d 41 m 10 s ) is located in the G333 
giant molecular cloud complex, at a distance of 3.6-kpc. 



It has a dust 
density of 2 x 



mass 
10 5 



of 



2.3 x 10 3 M w 
~ 3 dGarav et alj 



core was discovered by iGarav et al.l ( 2004 ) , comparing the 



and 
20041 ) 



an average 
The cold 



SES T Imaging Bolomet re Array (SIMBA) 1.2-mm sur- 
vey (jFaundez et al.ll2004o with Midcourse Space Experiment 
(MSX) mid-infrared data an d Infrared As tronomical Satel- 
lite (IRAS) far-infrared data. [Garay et all noted strong 1.2- 
mm continuum emission from this source in the absence of 
emission in any of the MSX and IRAS bands. From the 1.2- 
mm dust continuum, along with the upper limits of IRAS 
fluxes, they inferred th at the core is extremely cold (< 17 
K), massive and dense. IGarav et al.l hypothesised that the 
core is at an early stage of star formation, most likely before 
the development of an internal luminosity source, and will 
coll apse to form high-mass star(s). 

IPestalozzi et all (|2002h conducted SIMBA 1.2-mm ob- 
servations of the source based on the detection of a class II 
6.6-GHz methanol maser from survey by the Mt Pleasant 
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Observatory (|Ellingsen et al.l Il996l ). IPestalozzi et al.l sug- 
gested that there is a very deeply embedded object because 
of the detection of strong 1.2-mm emission and the methanol 
maser. Later, a class I 95.1-GHz methanol was detected by 
lEllingsenl (|2005h with the Mopra Telescope. 



2 OBSERVATIONS AND DATA PROCESSING 

The multi-molecular line data presented here were collected 
between July 2004 and October 2006 with the Mopra Tele- 
scope, operated by the Australia Telescope National Facil- 
ity (ATNF). It has a full width h alf-maximum (FW HM) 
beam size of ~32 arcsec at 100-GHz (|Ladd et al.ll2005l ). The 
observations were carried out with the new UNSW Mo- 
pra Spectrometer (MOPS) digital filterbank back-end and 
Monolithic Microwave Integrated Circuit (MMIC) receiver, 
except for 13 CO and C ls O which were observed with the Mo- 
pra Correlator (MPCOR) and the previo us Superconducto r 
Insulator Su perconducto r (SIS) receiver (|Bains et al.ll2006l ). 
According to lLadd et all (|2005l ) the main beam efficiency at 
86-GHz is 0.49, and at 115-GHz is 0.42. 

The observations are part of the on-going multi- 
molecular line survey (the 'Delta Quadrant Survey' or DQS) 
project carried out by the University of New South Wales, 
of the giant molecular cloud complex G333. The cold core 
G333. 125— 0.562 is part of the co mplex. The first pa per fea- 
turing 13 CO is already published (|Bains et al.ll200q ). In ad- 
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Figure 1. GLIMPSE 8-/im image (grey scale) overlaid with the 
SIMBA 1.2-mm dust continuum (red contours). The contour lev- 
els start from 690 mjy beam -1 with increments of 690 mjy 
beam -1 . The arrows show the position of the core and the IRAS 
source, the cross shows the 6.7-GHz methanol maser, while the 
scale bar denotes 2-pc at a distance of 3.6-kpc. No 8-^tm emission 
is evident at the position of the core. 
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Figure 2. The integrated emission map of SiO J = 2 — > 1 tran- 
sition (red contours) overlaid on CS J = 2 — > 1 transition (grey 
scale). The contour levels start from 2.4 K km -1 in steps of 0.8 
K km -1 (1 cr), in terms of the main beam brightness temper- 
ature, Tmb- The cross shows the 6.7-GHz methanol maser, the 
hatched circle indicates the beam size. The inset is the velocity 
profile of SiO averaged over the region indicated by the dashed 
box. Note the temperature scale of the SiO spectra is in antenna 
temperature. 



dition to the 3-mm multi-molecular line survey, NH3 and 
H2O (12-mm) maps of the core were also obtained during 
December 2006 with the Mopra Telescope. The main beam 
efficiency for 12-mm system is ??22GHz = 0.7 and FWHM 
beam size approximately 2 arcmin. The data presented in 
this work has a pointing accuracy of within ~ 5 arcsec. 



3 RESULTS AND DISCUSSIONS 

Shown in Figure Q] is the Spitzer Galactic Legacy Infrared 
Midplane Survey Extraordinaire (GLIMPSE) 8-fim image 
(grey scale) overlaid with the SIM BA 1.2-mm dust contin- 
uum (red contours) obtained from iMookeriea et al. I (|2004 ). 
From the image it can be seen that the core (bottom left) is 
isolated and has compact dust emission. No infrared emis- 
sion is evident from it. Inspection of the Spitzer MIPSGAL 
imagesat 24- and 70-pim show a source only at the longer 
wavelength. Thus this must be a cold dust core. From this 
core we have further detected emission from a number of 
molecular transitions, their observed parameters are sum- 
marised in Table [1] Among these, the detection of SiO (J 
= 2 — > 1) is of particular interest. Presented in Figure [5] 
is the integrated emission map of SiO (red contours) over- 
laid on CS (grey scale), both maps have been smoothed to 
36 arcsec. From the map we can see there are two strong 
peaks of SiO emission, one is the luminous IRAS source 
IRAS16172— 5028 and the other one is from the cold, dense 
core G333. 125—0.562. This is quite different to the extended 
distribution seen in other molecular transitions in this re- 
gion. Shown in the inset is the SiO velocity profile averaged 
over the core (dashed box), it is evident the line profile is 
broad. We will first discuss the distribution of the molecules, 
followed by the physical properties derived for the core. 



3.1 Distribution of the molecules 

Presented in Figure [3] are the integrated emission (zeroth 
moment) maps of the detected molecular lines over velocity 
range of -80 to -40 km s -1 , except 13 CO and C ls O, which 
are integrated over the ranges indicated in the figure. The 
molecular line emission peaks are coincident (within a beam) 
with the 1.2-mm dust, except for C ls O which is offset by ap- 
proximately 55 arcsec (~ 1-pc at a distance of 3.6-kpc). Note 
that there is extended C ls O emission towards the north of 
the dust peak. The low density tracer CO is clearly more 
extended than the other molecules, which ar e sensitive to 
denser gas. Recent studies (|Lintott et al.|[2005l ) suggest that 
CS is a good tracer of dense cores in high-mass star form- 
ing regions in contrast to low-mass star formation where CS 
may be depleted. Since the peak matches well with the dust 
peak, this suggests that CS does not suffer from significant 
depletion. However, we do note that depletion could hap- 
pen at smaller scales. Another feature evident in the HNC 
and HCO + maps is the tongue of extended emission to the 
south-west (arrow). A comparison between HCN and HNC 
integrated emission maps shows large differences in distri- 
bution, HCN being more extended than HNC, while at the 
same time the peak brightness of HCN is lower than HNC. 
Whether this is caused by optical depth effects needs further 
investigation. Although HCN and HNC have the same pre- 
cursors (HCNH + and H2CN + ) and similar dipole moments, 
it is believed that HCN is enhanced in warm environments 
in contrast to HNC, which forms preferentially in colder 
conditions (|Kim et al.ll2006l) . iHirota et all (|l998l ) found the 
abundance ratio of [HNC] /[HCN] rapidly drops as the tem- 
perature exceeds a temperature of 24-K. N2H + is known as a 
cold gas tracer as its abundance tends to be enhanced in cold 
dense cores, due to the depletio n of CO, which decreases the 
destr uction of N2H + and H^" (jWomack. Ziurvs fc Wvckofj 
ilmslevl 



ll992l ; lHotzel. Hariu fc Walmslevll2004l ). Rather than react 
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Figure 3. (a) to (j): Integrated emission maps of various molecular lines, as indicated for the core G333. 125—0.562. The contour levels 
start from 3 u of the r.m.s. The contour levels of N 2 H+, CH 3 OH, HNC, HCN, HCO+, C ls O and 1.2-mm continuum are in steps of 3 
a, and 1 a for CS, HC3N (both transitions), SO and SiO. The r.m.s. for CS and CH3OH is 1.1 K km s" 1 , N 2 H+ is 0.91, HNC, HCO+ 
and SiO are 0.82 K km s" 1 , HC3N (J = 10 - 9) and HCN are 0.61 K km s -1 , HC3N (J = 11 - 10) and SO is 0.91 K km s" 1 , and 
C ls O is 1.2 K km s" 1 . The temperatures are in TJj B . The integrated velocity range is from —80 to —40 km s — 1 , except C ls O which 
is from —80 to —54 km s — 1 . The arrows show the tongue extended emission in the south-west direction, (k): 13 CO integrated emission 
map from —62 to —54 km s , over a region covers the core (lower left) and IRAS16172-5028 (upper right). (1): Same as figure k but 
over velocity range —54 to —48 km s _1 . Contour levels for 13 CO maps start from 20 per cent of the peak, in increments of 10 per cent. 
The cross marks the peak of the 1.2-mm continuum and the square is the UCHII region associated with IRAS 16172—5028. The beam 
sizes are indicated with the hatched circle, the scale bar denotes 1-pc at distance of 3.6-kpc. 



ing with CO to form HCO + and H2, the rig*" ion reacts with 
N2 to form N2H + . From the integrated emission map of 
N2H + , we see a compact distribution of N2H + . This sug- 
gests the N2H + is tracing the outer cold envelope of the 
dense core. The other two molecules which show a com- 
pact distribution are CH3OH and HC3N (both transitions), 
both of which trace dense gas. SiO and SO are known 
to be greatly enhanced in outflows and shocked regions 
(jMartin-Pintado. Bachiller fc Fuentelll992t ). Detecting both 
these species suggests that energetic activity associated with 
shock waves is present. 



3.2 Physical properties 

In order to derive the physical parameters, spectra were spa- 
tially averaged over the FWHM angular size for molecules 
with compact emission, otherwise over a 48 arcsec region, 



then fitting with a Gaussian in CLAsfl The fitted parame- 
ters, i.e. line width, centre velocity, brightness temperature 
and optical depth (from hyperfine structure fitting of N2H + 
and NH3) were then used to derive column densities and 
masses. The derived parameters are summarised in Table 
[1] The column densities were calculated with an excitation 
temperature of 15-K, which is the average of the dust and 
NH3 values (see below). We note that the derived column 
densities vary by 20 per cent for the excitation temperature 
range (13 — 18 K). We have derived molecular abundances 
relative to the H2 column density, which is obtained from 
the 1.2-mm dust continuum, Nh 2 ~ 4 x 10 23 cm~ 2 . 

From the NH 3 (J,K) = (1,1) and (2,2) inversion 
transitions we have calculated the rotational tempera- 
ture (T12) of the core, following the met hod stated in 
lUngerechts. Winnewisser fc Walmslevl i|l986t ). We assumed 

Continuum and Line Analysis Single-dish Software, part of 
GILD AS software package by IRAM. 
jhttp://www.iram.fr/IRAMFR/GILDAS/^ 
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Table 1. Summary of the observed and derived parameters of the molecular lines: centre velocity ( V), line width (A V), peak brightness 
(Tpeak)i integrated brightness (J* Tmb dv), isotopic column density (AT), virial mass (M v ; r ) and FWHM angular size (9). The second 
last column is the abundance ratio (X) of the molecular line relative to H2. For lines with more than one velocity component, parameters 
reported here are for the component associated with the core velocity, —58 km s — 1 . The uncertainty range is in parentheses. 



Molecule 


Transition 


V 


A V 


-^"pcak 


J T MB dv 


N 


M viT 


8 


X 


Note 






(km s" 1 ) 


(km s _1 ) 


(K) 


(K km s -1 ) 


(10 14 cm" 2 ) 


(10 3 M ) 


(arcscc) 


(10- 10 ) 




3~ 111 111 






















13 co 


1-0 


-57.56 (0.03) 


4.90 (0.06) 


9.5 


49.3 (0.8) 


2.0 xlO 3 


1.2 


48 


5.0 x 10 4 


[ex] 


c 18 o 


1-0 


-56.63 (0.05) 


5.1 (0.1) 


2.6 


14.1 (0.3) 


5.2 xlO 3 


1.3 


18 


1.3 x 10 4 


[cx] 


cs 


2-1 


-58.8 (0.3) 


4.5 (0.7) 


0.9 


4.1 (1.3) 


1.7 


1.0 


18 


4.3 


[ex] 


c 34 s 


2-1 


-57.7 (0.3) 


3.1 (0.8) 


0.3 


1.0 (0.2) 


0.6 


0.5 


54 


1.5 


[c] 


HCO+ 


1 - 


-59.4 (0.1) 


4.6 (0.1) 


2.1 


10.4 (0.4) 


2.7 


1.1 


48 


6.8 


[ex] 


HCN 


1 — 






















F = 2 - 1 


-59.5 (0.1) 


2.2 (0.3) 


0.6 


1.4 (0.2) 














F = - 1 


-66.6 (0.1) 


5.3 (0.4) 


0.7 


3.9 (0.3) 














F = 1 - 1 


-51.8 (0.4) 
















[I] 


HNC 


1 - 


-59.07 (0.05) 


3.9 (0.1) 


2.2 


9.2 (0.5) 


> 3 x 10~ 3 


0.8 


46 


> 9 x 10~ 3 


[II] 


N 2 H+ 


1 - 


-58.8 (0.1) 


4.0 (0.1) 


1.0 


3.9 (0.3) 


24 


3.3 


99 


61 


[hf,c] 


CH3OH 


2 2 - loi A+ 


-58.0 (0.1) 


4.2 (0.2) 


1.2 


5.2 (0.3) 






84 




[c,III] 


HCCCN 


11 - 10 


-57.3 (0.2) 


3.6 (0.5) 


0.3 


1.3 (0.2) 


0.6 


1.3 


69 


1.5 


[c,IV] 


HCCCN 


10-9 


-57.41 (0.09) 


3.8 (0.2) 


0.6 


2.5 (0.1) 










[c] 


SO 


3 2 -2i 


-57.2 (0.2) 


5.0 (0.5) 


0.5 


2.7 (0.2) 


1.1 


1.3 


48 


2.6 


[11] 


SiO 


2-1 v=0 


-57.2 (0.3) 


9.0 (1.0) 


0.2 


1.8 (0.2) 


4.4 


3.9 


18 


11 


[II] 


12-mm 






















NH 3 


(1,1) 


-56.89 (0.02) 


3.46 (0.03) 


0.9 


10.7 (0.2) 


0.4 


1.5 


120 


1 


[hf,V] 


NH 3 


(2,2) 


-57.2 (0.1) 


3.7 (0.2) 


0.4 


2.7 (0.2) 










[hf] 


H 2 


616 — 523 


-53.9 (0.1) 


2.0 (0.3) 


0.1 


0.24 (0.03) 










[m] 



Note — 

[hf] stands for hyperfine structures; the peak brightness is for the main component, and the integrated brightness for the sum of 
all hyperfine components, except HCN, which we list each hyperfine component, 
[m] stands for maser. 

[ex] indicates the molecular line emission is extended, hence we used the angular size of dust continuum to derive the virial mass, 
[c] for compact emission. 

[el] indicates the emission is elongated, the angular size is the average of RA and DEC. 

[I] hyperfine component is blended with the —50 km s — 1 velocity component. 

[II] we derived the column density assuming the line is optically thin, hence it and the abundance ratio are lower limits. 

[III] no other CH3OH transition detected, therefore column density cannot be calculated. 

[IV] column density calculated using the rotational temperature (5.7-K) obtained from J = 10 — ► 9 and J = 11 — ► 10 lines. 

[V] the beam size at 22-GHz is ~ 2 arcmin, therefore the source is unresolved. We took the angular size to be one beam width as 
there is no noticeable extended emission. 



NH 3 (1,1) and (2,2) are tracing the same volume of gas, 
and all hyperfine components have the same excitation tem- 
perature. With the above assumptions, we have derived the 
rotational temperature, T12 = 13.1 K, with an unc e rtaint y 
range of 12.9 - 13.3 K. According to iDanbv et""aH ljl98&t ). 
below 20-K the rotational temperature of NH3 (l,l)-(2,2) 
matches the kinetic temperature (Tkin) well. 

We are also able to determine a dust temperature, mass 
and luminosity from the spectral energy distribution (SED). 
Inspection of the GLIMPSE and MIPSGAL images of the 
core only shows evidence of emission at the longest wave- 
length (70- /jm). Far-IR balloon-born e measuremen t s hav e 
been made at 150 and 210-^tm by iKarnik et all (|200ll ). 
through a 3 arcmin beam. Examination of the Spitzer in- 
frared images show this emission must be confined to the 
G333. 125— 0.562 core. Hence we use the 150-, 210-/xm and 
1.2-mm fluxes to determine the SED, constraining the size 
by the angular size (48 arcsec) of the 1.2-mm core, and ap- 
plying a greybod y fit with a dust emissivity index of (3 = 2 
(|Hill et al.ll2006l ). This yields the following parameters for 



the core: T dus t = 18.6 ± 0.1 K, L hol = (8.5 ± 0.3) x 10 3 
L and M gas = 1.8 x 10 3 Mq. Taken with the angular 
size these yield n H2 ~ 1.3 x 10 5 cm~ 3 and A r H 2 ~4x 10 23 
c m" 2 . We note these p arameters are similar to those derived 
bv lGarav et al.1 (|2004l ) who made use of just the 1.2-mm flux, 
combined with upper limits at 60 and 100-pm from IRAS 
data — yielding an upper limit to the temperature of 17- 
K. Thus, the dust-determined temperature is comparable to 
that derived from the NH 3 emission, giving us confidence in 
the conclusion that the core is exceedingly cold. If we as- 
sume that the molecular lines are thermalized then T^ n is 
approximately equal to the excitation temperature (Tex) for 
each molecule. 

Most of the detected molecular lines are peaked at ~ 
-57 km s-\ except for HCO+, HNC and N 2 H+ which are 
peaked at ~ —59 km s _1 . This likely due to lines being op- 
tically thick, for instance the H 13 CO + line peaks at —57 km 
s" 1 . The methanol maser also peaks at —57 km s -1 , with 
components at —53 and —63 km s _1 (|Ellingsenll20"05l ), in- 
dicative of outflows. The water maser is at —54 km s _1 . The 
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Figure 4. The integrated emission map of 13 CO (J = 1 — > 0) 
from —52 to —50 km s~ . The contour levels are 1.4, 4.3, 7.1, 
10, 21, 27 and 33 K km s^" 1 . The temperature is in terms of the 
main beam brightness temperature. The cross marks the cold core 
G333.125-0.562 and the square marks IRAS 16172-5028. The 
arrow indicates a likely region where the gas is being compressed. 



5 SUMMARY AND CONCLUSION 

From our molecular line survey of the G333 Giant Molecu- 
lar Cloud, we have detected thermal SiO emission from the 
massive, cold, dense core, G333. 125— 0.562, with an abun- 
dance enhanced over typical unshocked molecular cloud val- 
ues. This core has a gas mass of 1.8 x 10 3 Mq, is undetected 
up to 70-/im and has compact 1.2-mm dust continuum emis- 
sion. From the NH3 inversion lines we derived a temperature 
of 13 K, which is comparable with that derived from the SED 
(19 K). The detection of compact emission from a cold gas 
tracer (N2H + ) suggests this is a dense cold core. Typical 
line widths are between 4 to 5 km s _1 , indicating the core 
is turbulent, as expected in massive star formation. 

In conclusion, from these observations we believe that 
the cold massive core harbours a deeply embedded, massive 
protostellar object that is driving an outflow. This is oc- 
curring at a very early stage of star formation, prior to the 
creation of an infrared source in the core. 



line widths for most of the molecules are between 3 to 5 km 
s _1 as opposed to the ~ 0.5 km s _1 predicted for quiescent 
gas dominated by thermal broadening, suggesting the core 
is turbulent. The SiO however has the largest line width 
amongst the detected molecular lines, 9.0 ± 1.1 km s , 
leading further support to being associated with an outflow. 
Its abundance with respect to H2 has an upper limit of ~ 
10~ 12 in quiescent gas, but i s greatly enhanced in powerful 
shocks, up to ~ 10~ 7 (e.g. Ziurvs. Friberg fc Irvine! 1 19891 ; 



iMartin-Pintado. Bachiller fc Fuentelll992l ). We have deter 
mined the abundance to be ~ 10 -9 , in between these ex- 
tremes. A similar ab undance ratio has been found i n molec- 
ular outflows, e.g. ICodella. Bachiller fc Reipurthl (|l999h . 
iGarav et all (|2002T l. A follow up CO (J = 1 -> 0) molec- 
ular line observation over a five square arcmin region of 
the core was taken in June 2007; the spectra clearly shows 
line wings, however they are confused with multiple velocity 
components and need further investigation. 



4 WHAT IS GENERATING THE SIO 

EMISSION: OUTFLOWS OR COLLISIONS? 

The optically thick isotopomers clearly show two differ- 
ent velocity components, indicating there are two separate 
clouds in this region. The cloud at ~ —50 km s _1 is asso- 
ciated with IRAS 16172—5028, and is at the ambient ve- 
locity of G333 giant molecular cloud (see Figure |3j) . The 
other cloud at ~ —58 km s _1 is associated with the cold 
core G333. 125-0.562 (Figure EJc) . The SiO emission could 
be generated by these two clouds colliding, but we consider 
this unlikely. First, consider the confined emission of the 
SiO, localised to the dust core itself. If it is due to a cloud- 
cloud collision we would expect more extended distribution. 
Secondly, for a cloud-cloud collision there should be distinc- 
tive features in the density structure. Shown in Figure [3] is 
the 13 CO integrated emission map of velocity range —52 to 
—50 km s" 1 . Note the sharp edge (arrow) on the left side 
of the IRAS source (square) , suggesting compression of the 
gas here. However no compression is evident associated with 
the cold core (cross) at the velocity observed. 
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